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ABSTRACT 

The spatial distribution of compact dark matter in our Galaxy can 
be determined in a few years of monitoring Galactic globular clusters for 
microlensing. Globular clusters are the only dense fields of stars distributed 
throughout the three-dimensional halo and hence are uniquely suited to probe 
its structure. The microlensing optical depths towards different clusters have 
varying contributions from the thin disk, thick disk, bulge, and halo of the 
Galaxy. Although measuring individual optical depths to all the clusters is 
a daunting task, we show that interesting Galactic structure information can 
be extracted with as few as 40-120 events in total for the entire globular 
cluster system (observable with 2-5 years of monitoring). The globular cluster 
microlensing is particularly sensitive to the core radius of the halo mass 
distribution and to the scale length, surface mass density, and radial scale height 
variations of the thin disk. 

Subject headings: gravitational lensing, Galaxy: structure, globular 
clusters: general, Galaxy: fundamental parameters, Galaxy: halo, 
dark matter 



1. Introduction 

The search for compact dark matter in the mass range 10 _6 -10 2 M using microlensing 
(Paczyhski 1986) has come to fruition with the detection of many microlensing events 
towards the Galactic bulge, Large Magellanic Cloud, and Small Magellanic Cloud by various 
projects to monitor millions of stars (Alcock et al. 1993, 1996, 1997a,b; Aubourg et al. 
1993; Udalski et al. 1993). These studies have confirmed the presence of a massive bar in 
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the Galaxy by detecting more microlensing events than expected towards the Galactic bulge 
(Udalski et al. 1994; Alcock et al. 1995a; Paczyhski et al. 1994; Zhao, Spergel, & Rich 
1995) and have found evidence for a massive halo by detecting 8 times more events towards 
the Large Magellanic Cloud (LMC) than can be explained by known stellar populations 
(Alcock et al. 1997a). 

Uncertainties in the spatial distribution, kinematics, and masses of the lensing objects 
remain (cf. review by Paczyhski 1997). One of the major uncertainties is our incomplete 
knowledge of Galactic structure. Possible locations for the lenses include a maximal 
Galactic disk (Kuijken 1997), the LMC's disk (Sahu 1994), and the Galaxy's halo. Here 
we propose using globular cluster stars as background objects to discriminate the spatial 
distribution of the lensing objects statistically. This is possible because the globular clusters 
are distributed throughout the halo of our Galaxy, albeit concentrated towards the center. 
Therefore we can probe many lines of sight through varying amounts of the halo, disk, and 
bar of the Galaxy to determine the relative contributions of each to microlensing. Similar 
strategies using the ratios of Bulge, LMC and SMC optical depths have been suggested by 
Sackett and Gould (1993) to constrain halo flattening, and by Gould, Miralda-Escude & 
Bahcall (1994) to constrain relative contributions of the thin and thick disks and the halo. 

Another concern for microlensing studies is that the Galactic halo may be clumpy, 
as would be expected under the Searle and Zinn (1978) picture of Galaxy formation. 
The discoveries of the Sagittarius dwarf galaxy (Ibata, Gilmore & Irwin 1994) and of a 
possible intervening population towards the LMC (Zaritsky & Lin 1997) provide examples 
of clumps of stars, in tidal tails or dwarf spheroids in the halo, that might raise the number 
of microlensing events towards some lines of sight (Zhao 1997). Monitoring microlensing 
events towards many Galactic directions where globular clusters are located will average 
over such fluctuations in halo density and thereby alleviate the potential worry about 
non-representative lines of sight. 

Because the globular clusters are distributed throughout the Galaxy's halo, such a 
campaign of monitoring along different Galactic directions would be invaluable for weighing 
the compact massive component of the Galaxy (in stars or brown dwarfs) not only in the 
halo, but also in the thick disk, thin disk, and bulge. Measuring individual optical depths to 
the clusters is not practical, but interesting Galactic structure information can be extracted 
with as few as 40-120 events in total for the entire globular cluster system (observable with 
2-5 years of monitoring). 

In section 2, we review necessary background information on microlensing and on the 
globular clusters. We take M 15 as an example and discuss the observability of lensing 
events in globular clusters, using resolved stars and pixel lensing. In section 3, we study the 
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number of events required to distinguish between representative Galaxy models using Monte 
Carlo simulations. Finally, in section 4 we discuss some more implications of this work and 
summarize our conclusions. An appendix contains analytic results that supplement the 
simulations in section 3. 



The optical depth to microlensing is defined as the probability that a background 
source will be amplified by a factor 1.34 (Refsdal 1964; Vietri & Ostriker 1983). This 
amplification corresponds to the background source lying inside the Einstein radius of the 
lens. Because the area in the Einstein ring is proportional to lens mass, the optical depth 
depends only on the mean mass density and not on the masses of individual lensing objects, 
provided only that the Einstein radius substantially exceeds the physical size of both source 
and lens. This condition requires lensing masses M ^ 10 27 g ~ M e for length scales typical 
in the Galaxy. The geometry of source, lens, and observer affect the Einstein radius, and 
the full expression for the lensing optical depth r is (Vietri & Ostriker 1983) 



where D is the distance to the background source, L is the distance along the line of sight, 
p is the local mass density at distance L, and G and c are the gravitational constant and 
the speed of light. 

The Galaxy has ~ 150 globular clusters, distributed in a spheroidal halo about the 
Galactic center. Each typically contains ~ 10 4 -10 5 stars to a limiting magnitude of V < 22, 
giving a total of ~ 10 6 -10 7 stars to monitor. The optical depth to the clusters is 1CT 6 -1CT 8 , 
depending on the cluster location in the Galaxy and the Galactic model (Figure 1). In 
addition to the lensing optical depth, a full analysis would need to consider the mass 
distribution of the lenses and their velocities in order to predict the event rate and duration. 
These can then be compared to the time sampling of monitoring projects to determine 
the efficiency of a search. Such exhaustive analyses have been carried out for the Galactic 
bulge, Magellanic clouds, and M31 (Alcock et al 1995b), which are already being monitored. 
Observationally, the median duration of events towards the bulge is 30 days (Alcock et al 
1997a). The typical timescale for globular cluster microlensing events may be comparable 
or somewhat different, depending on the relative kinematics of globular clusters and the 
lensing populations. 



2. Microlensing of a Globular cluster 




(1) 



The number of stars in a globular cluster that can be resolved and individually 
monitored from the ground will depend upon the number of stars in the cluster, the degree 
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Fig. 1. — Histograms of microlensing optical depths to Galactic globular clusters under a 
"standard" model (model 5,6,7a in table 1). Left panel: Histogram of total optical depths. 
Right panel: Histogram of the optical depth due to the disk and bar components alone. The 
latter case is a lower bound to the expected optical depths, and omits the lensing contribution 
of both the stellar halo and the dark matter halo. 
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of concentration, the distance to the cluster and the seeing conditions. In the inner regions, 
crowding (up to 130 stars per square arcsecond for V < 22) is too severe for individual stars 
to be observed. Here one can instead monitor the brightness variations of pixels, which are 
a combination of many stars; this is the so-called pixel-lensing regime (Crotts 1992; Colley 
1995). 

We account for crowding by first calculating the effective number of stars that can be 
monitored in the cluster M 15, and then scaling the results for clusters of other luminosities 
and heliocentric distances. We chose M 15 as an example because stars in the interior core 
(down to 0.3") have been resolved by Hubble Space Telescope observations (Guhathakurta 
et al. 1996). M 15 is at a distance of ~ lOkpc, about the median distance to globular 
clusters. It has a high concentration index, which will lead us to a conservative (low) 
estimate for the fraction of resolved stars in a typical cluster. M 15 is also 2.7 times brighter 
(in V-band absolute flux) than the average cluster. 

The stars in the crowded interiors of Globular clusters can be monitored for pixel 
lensing events, where one pixel brightens due to the microlensing of one of the unresolved 
stars. The rate of detectable lensing events for each star is lower than in the uncrowded 
regime, because for a fixed amplification of the total flux in a pixel, a star needs to undergo 
a larger magnification, which depends on the fraction of light in the pixel arising from 
that star. This in turn depends on the number of stars in the pixel and the brightness of 
the star. Since the crowding varies continuously, decreasing from the unresolved core to 
resolved stars in the exterior, we calculate the mean effective optical depth for the entire 
cluster using a pixel microlensing formalism. In the outer regions of the cluster, where the 
surface density drops to £ 1 star per pixel, the lensing cross sections approach the typical 
Einstein ring area. 

Suppose a microlensing event is claimed for a fixed fractional increase in the brightness 
of a pixel. Following Colley (1995), 

Am = -2.5 x log ^ , (2) 

QJ* i Jpix 

where Am is the brightness change in magnitudes, A is the amplification, q is the fraction 
of the amplified star's light falling in the peak pixel of the point spread function, /* is the 
flux of the star, and f pix is the flux due to other sources in the pixel. The maximum radius 
u that produces amplification A is given by (Refsdal 1964) 

a= m ;+ 2 (3 ) 

where u = R/Re, and u 2 gives the area as a fraction of the area of the Einstein ring area 
Ti R\. If the star were resolved and the only one contributing to the flux in that pixel, 
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u — 1 and A = 3/v^5- Using the luminosity function of M 15 and its surface brightness 
profile (Guhathakurta et al 1996; Trager, King, Sz Djorgovski 1995) we calculate the loss of 
efficiency, which is the ratio of the area in the lens plane required to brighten the star in a 
field of surface brightness S vs when it is alone. This is a function of the surface brightness 
of the pixel and the brightness of the star. Summing over the stellar luminosity function 
f(m) for M 15 (Guhathakurta 1996) in V band, 



where u has an implicit dependence on the local surface brightness and hence on r. The 
value of £ varies from 6.4 x 10~ 4 at a radial distance of 0.1" from the center and approaches 
1 at r ~ 500" when stellar density reaches 0.12 stars per square arcsecond. With 3 x 10 5 
stars to a radius of 600" and a mean efficiency of {£) = 65.9%, the cluster can be treated 
like a cluster with 2 x 10 5 fully resolved stars. 

For other globular clusters, we retained the radial profile and luminosity function 
models from M 15. The total numbers of stars were scaled relative to M 15 according 
to each cluster's V band luminosity. The projected distribution on the sky (and hence 
crowding) and the number of stars above our V < 22 cutoff were adjusted for the distances 
to the individual clusters. The mean number of effectively resolved stars inside a 600" 
radius was 5.7 x 10 4 per cluster. 



We report here on a series of simulations to determine how well microlensing 
observations of globular clusters can distinguish among plausible Galactic models. In each 
test, the microlensing optical depth towards each cluster is calculated under two assumed 
Galactic models, and the probability of inferring the correct model is calculated for a given 
number of events. The results are then inverted to calculate the number of microlensing 
events needed to infer the correct Galactic model from the distribution of lensing events 
among the clusters with 90% confidence. 

The models generally consisted of three mass components: The thin disk, the halo, and 
the central bar. A fourth component, the thick disk, was added in one model. All the mass 
is assumed to be in compact objects. The optical depth to microlensing is calculated using 
equation (0), substituting the density p = p disk + p bar + p halo + p t hi c k disk- 

The models were required to obey three constraints. (1) Rotation Curve: The potential 
of the Galaxy consisting of these components was required to fit the rotation curve in the 




3. Distinguishing between Galactic models 
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inner Galaxy excluding the bar region (2.5 kpc < R < Rq) and a flat rotation curve with 
circular speed 220km/s in the outer Galaxy (R < R < 16.5 kpc) (Knapp, Tremaine, & 
Gunn 1978; Gunn, Knapp, & Tremaine 1979; Fich, Blitz, & Stark 1989; Malhotra 1994, 

1995) . Rq is the Sun-Galactic Center distance, and is taken to be 8.5 kpc throughout this 
work. (2) The second constraint is the local surface mass density of the thin disk, Eo, which 
lies in the interval 4OM /pc 2 £ S £ 8OM /pc 2 (Kuijken & Gilmore 1989, 1991; Bahcall, 
Flynn, & Gould 1992; Flynn & Fuchs 1994). We implement this by either fixing the local 
surface mass density or adjusting other parameters to ensure that the final value lies in the 
observationally acceptable range. This resolves some of the degeneracy between the radial 
force due to the disk and that due to the halo. (3) A third, implicit, constraint is the mass 
profile of the thin disk. From gas dynamics it is demonstrated to be an exponential with 
roughly the same scale length as the light (Knapp 1990; Malhotra 1995, 1994). 

The Thin Disk: The functional form of the disk density component (s) was taken to be 

p d = Eo/(2*i) x exp(-[r - Rol/n) x (sech^i]) 2 , (5) 

where T\ is the exponential scale length, z\ is the scale height, and (r, z) are Galactocentric 
cylindrical coordinates (Spitzer 1942). We took z 1 = 0.20 kpc at the Solar radius. The 
scale height was usually kept constant with radius. Such constancy is indicated by most 
studies of external galaxies (van der Kruit & Searle 1982), and is consistent with the 
diffuse near infrared light of the Milky Way (Spergel, Malhotra, & Blitz 1996; Freudenreich 

1996) . However, there may be exceptions to this rule (de Grijs & Peletier 1997), so in 
one case [model 5B] the scale height was allowed to increase with radius according to 
Z\ oc 1/S(r) oc exp(r/rx). This flaring is the appropriate behavior for a self-consistent 
isothermal thin disk with constant vertical velocity dispersion and exponentially decreasing 
surface density. 

Halo: A power law halo was adopted, following Evans (1993,1994) and Alcock et al 
(1995b): 

vl (l + 2g 2 )r c 2 + r 2 + (2-g-V 

Ph AnGq 2 (r 2 + r 2 + z 2/ q 2f ' ^ > 

where r c is the core radius of the halo and q the halo flattening parameter. Note that 
Alcock et al's parameter f3 has been set to (5 = 0, so that the asymptotic behavior at large 
radii is a flat rotation curve with circular speed vq. The power law halo component admits 
a self-consistent distribution function and an analytic potential; the latter was useful in 
ensuring that the models match observational rotation curve constraints. We do not include 
halo truncation, since relatively few globular clusters lie at galactocentric radii ^> 10 kpc 
while present evidence suggests the halo extends to 40-50 kpc (Little & Tremaine 1987). 
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Bar /Bulge: Detailed models of the central bar of our Galaxy have now been constructed 
by several groups (Dwek et al 1995; Zhao et al 1995; Stanek et al 1997). We chose to follow 
the Zhao et al (1995) bar model: 

n. - J^e"* 2 / 2 (7) 
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where a = 1.49 kpc, b = 0.58 kpc, and c = 0.40 kpc are the major axis, minor axis, and 
vertical scale lengths; s 4 = [{x/a) 2 + (y/b) 2 ] 2 + [zjcf\ and (x, y, z) is a Cartesian coordinate 
system aligned with the bar. We took the major axis of the bar to be at an angle with 
respect to the Sun-Galactic center line, and took the plane z = to be the Galactic plane 
(so that the bar is not tilted). In the rotation curve fits, we did not include data closer than 
2.5 kpc from the Galactic center, because of complex non-circular motions in the bar, and 
approximated the potential of the bar as a point source potential of the same mass. 

Tests: Seven tests were executed in all. Each determined the sensitivity of globular 
cluster microlensing observations to one parameter by comparing two paired models having 
a range of that parameter. Other structural parameters were adjusted as needed to maintain 
consistency with the rotation curve. The models are summarized in table 1. Test 1 is 
designed to discriminate between models with a large and small halo core radius. Other 
components necessarily become more massive in the case of a large core radius. Test 2 varies 
the halo flattening parameter, while keeping everything else fixed. Test 3 uses a fixed disk 
mass but a range of disk scale length; the halo parameters here change substantially also. 
Test 4 primarily explores local disk surface mass density. Test 5 allows the scale height of 
the thin disk to increase with Galactocentric radius. Test 6 includes a thick disk component 
of scale height 600 pc, surface density 1M Q / pc 2 , and scale length identical to the thin disk. 
Finally, test 7 varies the angle of the Galactic bar with respect to the Sun-Galactic center 
line. Note that the rotation curve constraint couples the parameters of different Galactic 
components in our models. Differences between microlensing optical depths in models A 
and B of a test may thus be due to changes in more than one Galactic component. 

For each model, the microlensing optical depth to each cluster was determined by 
evaluating equation ([]]) numerically for the density model given by equations (§)-(0). Our 
cluster sample consisted of 140 globular clusters whose coordinates, distances, and absolute 
V band magnitudes are tabulated in Djorgovski & Meylan (1993) and Djorgovski (1993). 
As a consistency check, we also implemented some of the halo models by Sackett & Gould 
(1993), and verified that our numerical integration code reproduced their optical depths 
towards the Magellanic Clouds and Galactic Bulge within 1%. 



Given the optical depths, we then approximated the effect of an observing campaign by 
defining the effective number n of trials per star monitored. This is essentially the number 
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of independent lensing events we expect to observe per unit optical depth, or equivalently, 
the number of event durations spanned by the observing campaign. Formally, n is defined 
by (E*) = nr, where (E*) is the number of events expected per star monitored and r is 
the lensing optical depth. Accounting for inefficiencies in the monitoring program, we write 
n m et/d. Here t is the duration of the monitoring program, d is the mean duration of 
lensing events, and e is the fraction of events that will be detected given the time sampling 
of the observing program. The observing efficiency e can be approximated as e = e^, 
where e± is the fraction of the year a source can be monitored and 62 is the fraction of 
events with durations longer than the time between observations and shorter than t. Note, 
though, that a precise calculation of n would require the full 6-dimensional phase space 
distribution function of the lenses to determine e. We use the approximations d ~ 30 days 
and e ~ 0.3, so that n ~ 3.5 x (t/yr). 

A simulation consists of assuming a uniform number of trials n for all clusters, and for 
each cluster drawing a random number of "observed" events from a Poisson distribution 
with mean (E) = N^nr under an assumed Galaxy model. The likelihood of obtaining the 
resulting fake data set is then computed under this model and under an alternate model. 
The model with the larger likelihood is taken as the inferred model. Many such simulations 
are run to obtain the probability of inferring the wrong model for each assumed model / 
alternate model pair. The test is then repeated for many values of n. In the limit of infinite 
data (n — > 00), the probability of error goes to zero. 

Table 1 contains values of the effective number of trials n$o required to distinguish 
between the two models in each test at the 90% confidence limit (i.e., with a 10% chance 
of inferring the incorrect model). Because this probability need not be the same for the 
two input models, two values of n 90 are tabulated for each test. To determine n 90 , we 
interpolate between measured error probabilities in simulations at a few values of n using 
scaling relations from appendix [A|: p(error) = erfc(/i)/2 where h oc yfn. These scaling 
relations can also be used to determine the values of n needed to achieve other confidence 
levels in distinguishing between models. 

These calculations do not include the optical depth due to self-lensing (i.e., lensing 
of one cluster star by another). This is justified by the low geometric weight given to the 
integrand in equation [L] when D — L ~ 1 pc. Dimensional analysis shows that for a cluster 
having mass M c \ and characteristic radius a, the self-lensing optical depth should be of order 
GM c i/(ac 2 ). We did more detailed calculations for clusters with density p oc (1 + (r/a) 2 ) 
for r < r t = af t and p = for r > r t . These yield mean optical depth to self-lensing 
(t s ) = GM d f(f t )/{2ac 2 ) = 2.4 x 10- 9 (M c/ /10 6 M Q )(l pc/a)/(r t ), where / « 3.5/ f t to 15% 
accuracy for 5 ^ f t <^ 500. For plausible cluster parameters, then, (r s ) <^ 10~ 9 , which is 
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negligible compared to foreground optical depths (Figure 1). 

The microlensing observations of the globular cluster system are most powerful in 
our tests 5 (flat vs. flaring disk), 1 (halo core radius), 3 (disk scale length), and 4 (local 
disk surface density). In general, the microlensing optical depths towards the globular 
clusters are sensitive to the large scale structural parameters of the Galaxy's most massive 
components, while they are not very sensitive to the presence or absence of a low mass 
thick disk (test 6) or to changes in the bar geometry within plausible limits (test 7). 
The sensitivity to test 5 is in part because model 5B (disk scale height increasing with 
galactocentric radius) is an extreme flaring model. We therefore take test 1 as our most 
favorable result for fully realistic models. 

Our simulations were constructed on the assumption that the same observational 
program is applied to each globular cluster. This assumption could be replaced by assigning 
different priorities to different clusters. The most reasonable observational program might 
leave out the poorest clusters entirely as providing an insufficient number of stars to be 
useful, while selectively including or rejecting other clusters depending on what facet 
of Galactic structure is being probed. For example, tests of disk flaring might omit 
observations of clusters at high Galactic latitude, where the disk contribution to r is 
negligible. Illustrative simulations with restricted subsets of clusters show that the duration 
of a monitoring campaign need only be doubled to achieve 90% confidence discrimination 
between model pairs with a well chosen sample of ~ 10 clusters. This strategy allows a 
specific model pair to be tested using a small fraction of the telescope time needed for a 
microlensing survey of the entire cluster system. Its results will, however, be more strongly 
affected by halo substructure. In addition, the optimum cluster subsample may be different 
for different model pairs. 



4. Conclusions 

From the simulations in this paper we see that the total number of events required to 
distinguish between models is typically ~ 40-120 for the whole set of clusters. 

The precise observing time required to achieve ~ 100 total events also depends on 
the lensing event duration. In this paper we have simply adopted a duration of 30 days, 
which is the median duration of microlensing events seen so far. The actual event durations 
will depend upon the masses and the spatial distribution of the lenses, and on the relative 
kinematics of the lens and the globular clusters in question. Statistics of event durations 
will, of course, provide further constraints on the masses, kinematics and spatial distribution 
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of the lenses. With this event duration of 30 days and monitoring the globular cluster 
system to a magnitude limit of V < 22 one should be able to distinguish between typical 
Galactic models discussed here [tests 1, 3, and 4] in 2-5 years. For comparison, ongoing 
microlensing projects have found 45 events in a year after monitoring ~ 10 7 stars towards 
the bulge and 8 events for ~ 10 7 stars towards LMC in 2 years. The event rate towards 
the globular cluster system should be between these two rates, because the clusters are 
distributed over a range of Galactic coordinates with a concentration towards the Galactic 
center. Two to five years is not an excessive amount of time given that the MACHO, 
OGLE, and EROS monitoring programs have been in operation for 5 years (1992-1997). 



Consider a comparison of two models under a fixed observing strategy. Let jjLj-i be the 
expected number of events and Tj :i the optical depth towards cluster j under model i. The 
two are related by = njN*jTj :i , where rij characterizes the observing time spent on 
cluster j (section |3|) and jV*^ is the effective number of resolved stars in cluster j (section |2|). 
Then if model 1 is correct, the probability of obtaining a data set whose likelihood under 
model 2 exceeds that under model 1 is 



A more readily calculated (though approximate) scaling is that P(error) « erfc (h) /2, where 
erfc(x) = 2 e~ x2 dxj^pn is the complementary error function, and where multiplying all 
the rij by a factor a changes h by a factor <Ja. Both these results are consistent with our 
simulations over the range of P (error) investigated. Interested readers may contact the 
authors for derivations of these results. 

We thank Jill Knapp, Ken Mighell, Bohdan Paczyhski and David Spergel for helpful 
discussions. JER's work is supported by a Kitt Peak Postdoctoral Fellowship. Kitt Peak 
National Observatory is part of the National Optical Astronomy Observatories, operated 
by the Association of Universities for Research in Astronomy. 



A. Analytic results for comparison of two models 
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Table 1: Parameters for tested Galaxy models. See text for definitions of model parameters. 
The units are as follows: S (M Q /pc 2 ); r\ (kpc); Vq (km/s); r c (kpc); M& (1O 1O M ); 
(degrees). Remaining quantities are dimensionless. The primary parameter tested in each 
model pair is in bold face in the table. The thin disk vertical scale height z\ = 0.20 kpc 
unless otherwise noted. lg[n 90 ] and E tot90 give the results of the simulations. E totj90 is 
essentially the number of lensing events required to distinguish a model from its alternative 
at the 90% confidence level; n 90 = E to t,9o/(J2 ^*,j T j) is the corresponding effective number 
of trials (as defined in text). The la random error in n 90 and E ta t,90 is 2.8%. 



*Disk scale height increases with Galactocentric radius. See text for description. 
**Thick disk present. Column density 1M Q / pc 2 at Sun; scale height 0.60 kpc, scale length 3 kpc. 



